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Abstract 

We present a measurement of CP-violating asymmetries in B° — > it + it~ decays based on a 
41.8 fb -1 data sample collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. We fully reconstruct one neutral B meson as a B° — > tt + it~ 
CP eigenstate and identify the flavor of the accompanying B meson from its decay products. 
From the asymmetry in the distribution of the time intervals between the two B meson decay 
points, we obtain the CP-violating asymmetry parameters S n7T = —1.21 ^^(stat) iai3( s y s t) an d 
Aw = +0.94 1^31 (stat) ± 0.09(syst). 

PACS numbers: PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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Kobayashi and Maskawa (KM) proposed, in 1973, a model where CP violation is incor- 
porated as an irreducible complex phase in the weak- interaction quark mixing matrix [[!]. 
Recent measurements of the CP-violating parameter sin20i by the Belle and BaBar || 
collaborations established CP violation in the neutral B meson system that is consistent 
with KM expectations. Measurements of other CP-violating parameters provide impor- 
tant tests of the KM model. In this Letter we describe a measurement of CP-violating 
asymmetries in the mode B° — > 7r + 7r~ ||; these are sensitive to the parameter sin20 2 ||. 

The KM model predicts CP- violating asymmetries in the time-dependent rates for B° 
and B° decays to a common CP eigenstate, fcp [Bfl- In the decay chain T(4S) — ► B°B° — > 
fcpftag, where one of B mesons decays at time tcp to fcp and the other decays at time 
t ta g to a final state /tag that distinguishes between B° and B°, the decay rate has a time 
dependence given by |7j 

e -|At|/T B o 

P^(At) = — [1 + q ■ {S^ sm(Am d At) 

At b o 

+A 7T7T cos(Am d At)}} , (1) 

where r^o is the B° lifetime, Am d is the mass difference between the two B° mass eigenstates, 
At = tcp — ttag, and the 6-flavor charge q = +1 (—1) when the tagging B meson is a 
B° (B°). The CP- violating parameters «S 7r7r and A nn defined in Eq. (|T|) are expressed by 
S n7T = 2PmA/(|A| 2 + 1) and A^ = (|A| 2 — 1)/(|A| 2 + 1), where A is a complex parameter that 
depends on both B°B° mixing and on the amplitudes for B° and B° decay to 7r + 7r~. In the 
SM, to a good approximation, |A| is equal to the absolute value of the ratio of the B° to B° 
decay amplitudes. We would have S nn = sin20 2 and A^ = 0, or equivalently |A| = 1, if 
the b — ► u tree amplitude were dominant. The situation is complicated by the possibility of 
significant contributions from gluonic b —>■ d penguin amplitudes that have a different weak 
phase and additional strong phases ||. As a result, S nn may not be equal to sin 202 and 
direct CP violation, A nn ^ 0, may occur. 

This measurement is based on a 41.8 fb _1 data sample, which contains 44.8 million 
BB pairs, collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider operating at the T(45) resonance. At KEKB, the T(4S I ) is produced 
with a Lorentz boost of /3^y = 0.425 nearly along the electron beamline (z). Since the B° 
and B° mesons are approximately at rest in the T(4S) center-of-mass system (cms), At 
can be determined from the displacement in z between the fcp and / tag decay vertices: 
At ~ (z C p - z tag )/f3>yc = Az/(3^c. 

The Belle detector [TI| is a large-solid-angle spectrometer that consists of a silicon vertex 



detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a super- conducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented 
to detect K° L mesons and to identify muons (KLM). 



The B° — ► 7r + 7r~ event selection is described in detail elsewhere ||11|| . We use oppo- 
sitely charged track pairs that are positively identified as pions according to the combined 
information from the ACC and the CDC dE/dx measurement. Candidate B mesons are re- 
constructed using the energy difference AP = Eg™— ^beam anc ^ the beam-energy constrained 



mass Mb c = y(Pbc^ s m ) 2 — (Pb 118 ) 2 , where Pbeam i s the cms beam energy, and Pg ms and p^" s 
are the cms energy and momentum of the B candidate. The signal region is defined as 5.271 
< Mb c < 5.287 GeV/c 2 and |AP| < 0.067 GeV, corresponding to ±3cx from the central 
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values. In order to suppress background from the e + e — > qq continuum (q = u, d, s, c), 
we form signal and background likelihood functions, Cs and Cbg, from two variables. One 



is a Fisher discriminant determined from six modified Fox- Wolfram moments [|12|; the other 



is the B flight direction in the cms, with respect to the z axis (cos#b). We determine 
Cg from Monte Carlo (MC) and Cbg from data, and require £s/(£s + £-bg) > 0.825 for 
candidate events. Figure [I] shows the AE distribution for 7r + 7r~ candidates. The signal 
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FIG. 1: AE distribution for tt + tt event candidates that are in the M^ c signal region. 

yield is extracted by fitting the AE distribution with a Gaussian 7r + 7r~ signal function, plus 
contributions from misidentified B° K + tt~ events, three-body i?-decays, and continuum 
background. From the fit, we obtain 73.5 ± 13.8(stat) 7r + 7r~ events, 28.4 ± 12.5(stat) K + tt~ 
events, and 98.7 ± 7.0(stat) continuum events in the signal region. The K + n~ contami- 
nation is consistent with the K — > n misidentification probability measured independently. 
The contribution from three-body 5-decays is negligibly small in the signal region. 

Leptons, charged pions, and kaons that are not associated with the reconstructed B° — > 
7r + 7r~ decay are used to identify the flavor of the accompanying B meson. We use two 
parameters, q and r, to represent the tagging information. The first, q, is already defined 
in Eq. (|1|). The parameter r is an event-by-event, MC-determined flavor-tagging dilution 
factor that ranges from r = for no flavor discrimination to r = 1 for unambiguous flavor 
assignment. It is used only to sort data into six r intervals (boundaries at 0.25, 0.5, 0.625, 
0.75 and 0.875). The wrong tag fractions for the six r intervals, wi (I = 1,6), are determined 
from the data; we use the same wi values that were used for the sin20i measurement 0. 

The vertex positions for the ii + 7i~ and /tag decays are reconstructed using tracks with 
associated hits in the SVD. Each vertex position is also constrained by the interaction point 
profile smeared in the r-0 plane by the average transverse B meson decay length. The / tag 
vertex is determined from all well reconstructed tracks, excluding the B° — > 7r + 7r~ candidate 
and tracks that form a Kg candidate. 

The proper-time interval resolution for the signal, R sig (At), is obtained by convolving a 
sum of two Gaussians (a main component, plus a tail component caused by poorly recon- 
structed tracks) with a function that takes into account the cms motion of the B mesons. 
We use the same parameters as those used for the sin20i measurement; the fraction and 
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the typical width of the main Gaussian are 0.97 ± 0.02 and 1.49 ps, respectively [Q]. The 
background resolution function R q g(At), which is dominated by continuum background, has 
the same functional form but the parameters are obtained from a sideband region in Mb c and 
AE. Using these resolution functions, we perform a B° lifetime measurement that yields 
T£o = 1.49 ± 0.21 (stat) ps for B° — > 7r + 7r~ candidates, which is consistent with the world 
average value |T3] ]. 

We determine CP violation parameters by performing an unbinned maximum-likelihood 
fit of a CP-violating probability density function (pdf) to the At distributions. We define 
the likelihood value for each event as a function of «S W7r and A^: 

Pi = J \{f^VUAt\ wv, S™, A nw ) + f l Kw V q K7T (At>, wi)} ■ R sig (AU - At') 

+f l q gV q - q (At') ■ R m {AU - At')\dAt'. (2) 

Here f^, f l Kn , and f qq (= 1 — f l n7T — f l K7T ) are the fractions of n + 7r~ signal, K + n~ background, 
and continuum background in flavor-tagging interval /, respectively. These fractions are 
determined on an event-by-event basis as a function of AE and Mb c , properly normalized 
by the average signal and background fractions in the signal region. The average fractions 
of qq background for six r bins (I = 1, 6) are 0.632, 0.505, 0.462, 0.440, 0.322 and 0.117. For 
higher r values where we are more sensitive to the asymmetry, the fraction of continuum 
background decreases; the ratio of 7r + 7r~ signal events to background K + ir~ events is the 
same for all r bins. The pdfs for 7r + 7r~ (V% w ), K + n~ (V 1 ^^, and continuum background (V q q), 
are convolved with their respective resolution functions. We use the same vertex resolution 
function for 7r + 7r~ and K + it~ candidates. For the 7r + 7r~ signal, the pdf is given by Eq. ([!]) 
with q replaced by q(l—2wi), to account for the dilution due to wrong flavor tagging. The pdf 
for the K+-K- background is V^(At, w t ) = e- |A * l/r s°/4r B o{l+g-(l-2w i )^l A - 7r cos(Am d At)}, 
where Ak-k is the B° — > K~n + and B° — > K + ir~ decay rate asymmetry. We fix Ak-k — fllT] , 
and tbo and Am^ to their world average values jl3| . The pdf used for the qq background is 
V q q(At) = {/ T e~ |A * l/Tbkg /2r b kg + (1 - f T )6(At)}/2, where f T is the background fraction with 
an effective lifetime r^ g and 5 is the Dirac delta function. We determine f T = 0.011 ± 0.004 
and Tbkg = 2.7io;? ps from the sideband data. In the fit, S W7r and A T7T are free parameters 
determined by maximizing the likelihood function £ = Hi Pi, where the product is over all 
B° — > 7t + tt~ candidates. 

The result of the fit to the 162 candidates (92 B°- and 70 5°-tags) that remain after 
flavor tagging and vertex reconstruction is: 

S n7T = -1.21 t^?(stat) lg;^(syst); 
A** = +0.94 1^ (stat) ±0.09(syst). 

The result is 1.3<r away from the physical boundary + A\ v = 1, which is consistent 
with a statistical fluctuation. The correlation between S nn and A^n is found to be 0.28. In 
Figs. 0(a) and (b), we show the At distributions for B°- and £?°-tagged events together with 
the fit curves; the background-subtracted At distributions are shown in Fig. 0(c). Figure 0(d) 
shows the background-subtracted CP asymmetry between the B°- and _B°-tagged events as 
a function of At. The result of the fit is superimposed and shown by the solid curve. 

The systematic error on S nn is primarily due to uncertainties in the background fractions 
(±0.09) and a possible fit bias near the physical boundary (to.H)- F° r A n7T , the background 
fractions (±0.06) and the wrong-tag fractions (±0.06) are the two leading components. We 
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FIG. 2: The At distributions for the B° — > 7r + 7r~ candidates in the signal region: (a) candidates 
with q = +1, i.e. the tag side is identified as -B ; (b) candidates with q = — 1; (c) ir + ir^ yields 
after background subtraction. The errors are statistical only and do not include the error of the 
subtracted background obtained by a fit. The rightmost (leftmost) bin ranges from 5 to 10 ps (—5 
to —10 ps); (d) the CP asymmetry for B° — > tt + tt~ after background subtraction. The point in 
the rightmost bin has a large negative value that is outside of the range of the histogram; (e) the 
raw asymmetry for B° — > tt + tt~ sideband events. In Figs, (a) through (c), the curves show the 
results of the unbinned maximum likelihood fit. In Fig. (d), the solid curve shows the resultant 
CP asymmetry, while the dashed (dotted) curve is the contribution from the cosine (sine) term. 



find that uncertainties in f^ w , f l Kn and f qq in Eq. (0) account for the largest systematic 
errors. We add their contributions in quadrature to obtain the above values, where each 
contribution is obtained by varying a parameter by its error and repeating the fit. Other 
sources of systematic error are uncertainties in the resolution function, physics parameters 
(Amrf, t b o, and Ak-k) and the background modeling. A value of Ak-k = —0.06 ± 0.08 is 
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obtained from the self-tagged B° — > K + tt sample [|TTJ; this introduces a systematic error 
of < 0.01 for Smr and igio? for Aw- 

We validate our fitting procedure using a large ensemble of MC pseudo-experiments 
wherein events are generated with nominal pdfs and the observed number of events. For 
various input values of S nn and A nn , we confirm that there is no bias in the fit results. 
We find that the average expected errors, 0.52 for S W7r and 0.35 for A nn , are larger than 
our measurements. However, the probability of obtaining errors that are smaller than our 
measurement is 5.4% for S nn and 23.6% for A^\ the results are within the expected range 
of statistical fluctuations. The measured likelihood is in good agreement with the average 
likelihood value obtained in the pseudo-experiments. 

We perform a number of cross checks. We examine the event yields and At distributions 
for B - and 5°-tagged events in the sideband region and find no significant asymmetry as 
shown in Fig. ^|(e). We select B° — > K + 7t~ candidates, which have the same track topology 
as B° — > 7r + 7r~, by positively identifying charged kaons. A fit to 309 candidates (209 ± 16 
signal events) yields t b o = 1.73 ± 0.15(stat) ps and Amj = 0.57 ± 0.08(stat) ps _1 ; these 
are consistent with the world average values |R|. Ak-k is 0.07 ± 0.17, in agreement with 
the counting analysis mentioned above and Sk-k — 0.15 ± 0.24, which is consistent with 
zero. We also select B° — * D~tt + , D*~tt + and D~p + candidates using the same event 
shape criteria. Neither mixing-induced nor direct CP- violating asymmetry is observed as 
expected. As an additional test of the consistency of the background treatment, we add 
events from the B° — > 7r + 7r~ sideband and adjust their AE and M^ c values. A fit to this 
background-enriched control sample, which has a background fraction comparable to the 
B° — > 7r + 7r~ sample, yields S = 0.08 ±0.06 and A = 0.03 ±0.04, both consistent with a null 
asymmetry. 

We determine the statistical significance from the likelihood function, taking into account 
the boundary of the physical region as well as the effect of the systematic error. The Feldman- 
Cousins frequentist approach ]TJ| gives a 99.6% confidence level (C.L.) for —1 < < 0, 
equivalent to a 2.9cr significance for a Gaussian error. A similar analysis yields a significance 
of 2.9a for 0<i„< 1. The 95% C.L. intervals are found to be -1.00 < S n7T < -0.39 and 



±0.30 < A™ < ±1-00, respectively [15 



In summary, we have measured the CP violation parameters in B° — > 7r + 7r _ decay. Our 
result for S n7T indicates that mixing-induced CP violation is large. The large A V7T term is 
an indication of direct CP violation in B meson decay, and suggests that there is a large 
hadronic phase and interference between the tree and penguin amplitudes fllfj. In this case 



the precise determination of sin 202 from «S T7r requires additional measurements including 
the branching fractions for the decays B° — > 7r°7r° [17 . 
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